How viscous is a superfluid neutron star core? 



N. Andersson*, G. L. Corner^ and K. Glampedakis* 
* School of Mathematics, University of Southampton, Southampton SOU IB J, UK 
''Department of Physics, Saint Louis University, St. Louis MO 63156-0907, USA 

(Dated: February 2, 2008) 

We discuss the effects of superfluidity on the shear viscosity in a neutron star core. Our study 
combines existing theoretical results for the viscosity coefficients with data for the various superfluid 
energy gaps into a consistent description. In particular, we provide a simple model for the electron 
viscosity which is relevant both when the protons form a normal fluid and when they become 
superconducting. This model explains in a clear way why proton superconductivity leads to a 
significant strengthening of the shear viscosity. We present our results in a form which permits the 
use of data for any given modern equation of state (our final formulas are explicitly dependent on 
, the proton fraction). We discuss a simple description of the relevant superfluid pairing gaps, and 

construct a number of models (spanning the range of current uncertainty) which are then used to 
' discuss the superfluid suppression of shear viscosity. We conclude by a discussion of a number of 

challenges that must be met if we are to make further progress in this area of research. 



I. INTRODUCTION 



Neutron stars arc often thought of as exciting cosmic laboratories. This is natural since their description depends 
on much complex physics [| : With a mass of about one and a half times that of the Sun compressed inside a radius 
, of ten kilometers or so, they are compact enough to require a fully general relativistic description. With central 
■ densities several times the nuclear saturation density it may be energetically favourable for exotic phases of matter, 
OO \ like kaon condensates, hyperons and/or deconfmed quarks, to be present in the core. With temperatures much below 
the Fermi temperature for the various constituents, neutron stars are cold on the nuclear scale which means that 
the presence of both solid (the outer layers form a kilometer sized nuclear lattice) and superfluid regions is expected, 
j . Neutron star dynamics is also intriguing. The obvious example of this is the glitches, which are taken as evidence of 
at least two weakly coupled interior components, observed in a number of radio pulsars Q- In addition, the evolution 
of the sample of the neutron star population which is spinning rapidly may be affected by both hydrodynamical 
and/or radiation-driven instabilities [3j. This is particularly interesting from the point of view of gravitational- wave 
observations. With a generation of large-scale gravitational-wave interferometers now reaching design sensitivity, it 
I ■ is worth emphasising that a signal from an (unstable?) neutron star pulsation mode could provide a unique, and 
perhaps the only, probe of the bulk motion of matter inside such stars. The information that could be gleaned from 
,^3 . such data would be truly unprecedented 0, |f| . 

As a useful illustration of a problem which plays a key role in neutron star dynamics, let us consider the coupling 
between the elastic crust of nuclei and the fluid core. The relevant coupling timescale impacts on, for example, the 
damping rate of oscillations driven unstable by gravitational-radiation emission (eg . the r- modes (3), the relaxation 
following a pulsar glitch |6|, LJ |8( , and possible neutron star free precession |9|, llOt lll| . In its simplest form the problem 
is analogous to the classic "spin-up problem" in fluid dynamics, where one studies the rotation rate of a viscous fluid 
following a change in rotation of the container. It is well-known that the fluid velocity changes due to the formation 
of a so-called Ekman layer at the interface. The role of this viscous boundary layer is to ensure that the fluid motion 
satisfies a no-slip condition at the interface. Its presence results in a coupling to the bulk of the fluid on a timescale 



< B ~*,(|) with (1) 

where 5e is the width of the Ekman layer [g^, R is the radius of the container, Q is the rotation rate, p is the mass 
density at the interface and 77 is the shear viscosity coefficient. For a neutron star, where 5e/R ~ 1CP 6 it is clear that 
the resulting timescale is much shorter than that of viscous diffusion, i.e. 

(2) 

V 

Analogously, in the case of an oscillation in the fluid the Ekman layer provides an efficient damping mechanism, see 
Q for a recent review. 

However, it is a serious oversimplification to model the core-crust interface as a viscous fluid/solid wall transition. 
In particular, the fact that the crust is expected to be permeated by superfluid neutrons should also be considered. 
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Similarly, the multi- fluid nature of the expected superfluid neutron/superconducting proton mixture in the outer core 
must be accounted for. This involves many difficult issues associated with the presence of rotational vortices in the 
neutron superfluid, and magnetic fluxtubes in the proton fluid (if it forms a tv pe-II superconductor). The interactions 
involving the vortices, leading to the non-dissipative entrainment effect |l2j as well as the mutual friction |l3l |. and 
the potential pinning of vortices to the nuclear lattice [7], |8( are not yet understood in detail. It is clear that one faces 
a problem which stretches our understanding of neutron star physics considerably. In a series of papers, we aim to 
make progress by considering each of the relevant issues in turn. 

The purpose of the present paper is to discuss the shear viscosity in a superfluid neutron star. In doing this, we 
are taking a small step towards a fully consistent description of dissipative neutron star cores. As a first step in this 
development we revisit the standard shear viscosity (ignoring for the moment the multi- fluid dynamics aspects), and 
ask how well understood the relevant viscosity coefficients are. This leads to a set of useful formulas that can readily 
be used in future work. In addition, our discussion identifies a number of challenging research problems that require 
further attention. 



II. REVISITING THE SHEAR VISCOSITY 



We begin by making the following observation: In their by now classic paper on the effect of viscosity on damping of 
oscillations in superfluid neutron stars, Cutler and Lindblom |l4j state: "It is interesting to note that contrary to our 
experience with other superfluids like He4, neutron star matter becomes more viscous in the superfluid state than it 
was in the normal state." This, at first sight, counter-intuitive result provides ample motivation for our investigation. 

Estimates of the effect that shear viscosity has on the fluid motion in neutron stars have so far almost exclusively 
been based on the work of Flowers and Itoh |15llT^| . who calculated the relevant transport coefficients at supranuclear 
densities. Based on the results of these papers, Cutler and Lindblom 0] discussed the viscous damping of oscillations 
in the neutron star fluid. Their analysis was based on, first of all, a fit to the result in ^(| for the total shear viscosity 
in normal fluid neutron stars, which can be written 

Vn a 2 x IQ 20 P % A T- 2 g/cm s (3) 

where P15 = p/10 15 g/cm 3 and Tg = T/10 8 K. As stated by Flowers and Itoh (and as is also evident from Figure 2 
in |16|). the dominant contribution to r\ n is due to neutron- neutron scattering. This result makes sense since the 
neutrons make up the bulk of the fluid at the relevant densities (and more exotic massive particles like hyperons are 
not being considered) . However, the neutrons in the outer core are expected to become superfluid as soon as the star 
cools below (say) 10 9 K, i.e. soon after its birth in a supernova core collapse. Below the neutron superfluid transition 
temperature the neutron-neutron scattering is suppressed and the dominant contribution to the shear viscosity is 
made by scattering processes involving the relativistic electrons. Cutler and Lindblom 0] argue that the relevant 
estimate (for electron-electron scattering) in a superfluid neutron star is 

r/ e c e L « 6 x 10 2( V 15 r 8 - 2 g/cm s (4) 

This estimate for r) ee has since been used in a number of contexts. The puzzling fact that rj^ > rj n throughout a 
typical neutron star does not seem to have caused much concern. In addition to possible confusion about this result, 
the estimate (0J is not quite satisfactory. One reason is that it is associated with a particular supranuclear equation 
of state (derived by Baym, Bethe and Pethick in 1971 0]). I* 1 practise, one would like to have an estimate which 
explicity includes the dependency on the electron (=proton) fraction. This is important since one may then be able 
to better understand the effect of varying the supranuclear equation of state (and the proton fraction) within the 
current range of uncertainty. 

To derive the relevant estimate, we take as our starting point the general formula (cf. 0,0,0]) 

n ePl , K s 

where n e is the electron number density, p e = Tikp — ?i(37r 2 n e ) 1 / 3 is the corresponding Fermi momentum, m* — p e /c 
is the effective electron mass and r is the timescale for momentum transfer due to the various scattering processes in 
which the electrons are involved. In a normal fluid neutron star core, scattering off of the (nonrelativistic) protons 
provides the main channel for momentum exchange. However, if the protons become superconducting then electron- 
electron scattering takes over as the dominant contribution to 



it the pr< 



3 



In order to model the resultant shear viscosity, we use the timescale for electron-proton scattering derived by Easson 
and Pethick 111: 



' ep 99 
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(6) 



where a = 1/137 is the fine-structure constant, tp p is the proton Fermi energy and ku is the Thomas-Fermi screening 
wave- vector (to be discussed later). Note that we are assuming that the matter is charge neutral, i.e. we take n e = n p . 

We can extract an analogous formula for electron-electron scattering from Eq. (91) in This equation can be 
written 
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Noting that 
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where x p is the proton (electron) fraction, we see that we only need to retain the leading order contribution (in 
m e c/p e ). By comparing (JSj and J7J we then find that 
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where CFe — hckp for relativistic electrons. 

In order to combine the two results for r ep and r ee , we need to note that individual scattering processes add like 
"parallel resistors" , cf. , which means that we have 
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This shows that the most important contribution to the shear viscosity comes from the most frequent scattering 
process (as one would intuitively expect). Using 
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where €f p — h 2 k F /2m* and m* is the effective proton mass, we find that r ee >> r ep under normal circumstances. 
Thus the electron-proton scattering provides the dominant contribution to the electron shear viscosity, in accordance 
with the discussion of Flowers and Itoh 0] . 

We now want to account for the likely possibility that the protons are superconducting. In that case the electron- 
proton scattering will be suppressed, essentially because there will be fewer states available for the protons to scatter 
into. In order to allow for the transition to proton superconductivity, we introduce a suppression factor 1Z P such that 



' ep 
Tlr, 
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Far below the critical transition temperature at which the protons become superconducting we should have TZ P 
and we see from (11011 that the electron-electron scattering then dominates the shear viscosity @. 
Combining the above results we obtain the following final formula 
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(13) 



which should be used together with It should be noted that we have neglected the contribution from electron- 
neutron scattering. This should always be a valid approximation. 
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Let us now turn to the screening factor kjt. Electron scattering is screened by both electrons and protons, and the 
relevant screening wave vector is given by [66j 

The first term in the bracket is due to the electron screening, while the second is due to the protons. It is straight- 
forward to show that the latter tends to dominate under the conditions that prevail in a neutron star core. 

Having introduced the key ingredients, we can discuss two limiting cases. In the case of normal protons we take 
n T = 1 in fT3|l. This leads to 

Vep « 1.8 x 10 18 (^_) pTt 8 - 2 g/cm s (15) 
Here, and in the following, we have taken (for simplicity) m* = m p in evaluating kf t . Our estimate l|15|) agrees with 

< 03. 



the result derived by Easson and Pethick |18|. In the opposite limit the protons are (strongly) superconducting, which 
means that 1Z P — 0, and we arrive at the estimate 

Vee « 4.4 x 10 19 ( J581) P %%-* g/cm s (16) 

This result should be compared to the formula used by Cutler and Lindblom |14| in deriving J3J . One can easily show 
that the two results are compatible. This means that the differences that can be seen in Figure ^ where we compare 
the various coefficients, is entirely due to the fact that different equations of state (eg. proton fractions) are used. 
This figure shows clearly that the (potentially puzzling) result that the electron-electron viscosity in the neutron star 
core is stronger than that due to normal neutron-neutron scattering remains true. We also learn that the key reason 
for the emergence of the dominant electron shear viscosity is the superconductivity of the protons. 

The results illustrated in Figure ^ (and others discussed throughout the paper) were obtained using one of the 
phenomenological PAL equations of state (with compression modulus K — 240 MeV) .22]. In order to obtain a 
suitably simple model we have followed Kaminker et al \2?\ who provide the following fit for the two nucleon number 
densities; 

n x = ap\j (1 + cpi 4 + dpj 4 ) (17) 

where 



a = 0.1675, 6=1.8185, c = 2.0288 , d = 0.02444 for x = n 

and 

a = 0.0006823, 6 = 2.6727, c = 0.1946 , d = 0.01604 for x = p 

The total number density rij, and the proton fraction x p follow immediately from n.& = n n + n p and x p = n p /rib, 
respectively. 

In one of the few alternative calculations of the neutron shear viscosity (that we are aware of) , Mornas [24j derives 
the following expression [cf. her Eq. (4.9b)] 

r, n « 2.4 x 10 19 pi^T 8 - 2 g/cm s (18) 

It is notable that the scaling with the density is very different from that in Eq. J3J), possibly because of the density 
dependence of the effective neutron mass. As can be seen from Figure ^ Eq. I|18|) leads to a considerably weaker 
viscosity at high density. It should also be noted that Mornas's formula leads to results that are very close to the 
electron-electron scattering results for the crust region determined by Flowers and Itoh 15]. We have no rational 
explanation for why this should be the case (in fact, it must be a coincidence), but it is a potentially useful observation. 



III. MODELLING THE EFFECTS OF SUPERFLUIDITY 



The formulas given in the previous section provide a useful step towards modelling the dynamics of realistic neutron 
stars since they allow us to determine a consistent electron viscosity for various supranuclear equations of state and 
core temperatures. Of course, in order to use these formulas we need not only the total number density and the proton 
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FIG. 1: An illustration of the various shear viscosity estimates for a temperature of 10 s K. The thick solid line shows our 
estimate -q ee while the thick dashed curve corresponds to rj ep (where we have taken mp 1 = m p in evaluating the screening 
factor). These estimates are compared to i) the result of Cutler and Lindblom [Til ] for electron-electron viscosity (CL87, thin 
solid line), ii) the total shear viscosity in a normal fluid neutron star, which is dominated by neutron-neutron scattering |T3. fl^| 
(FI79/CL87, thin dashed line), iii) an alternative calculation of the neutron shear viscosity due to Mornas [24| (M94, dot-dashed 
thin line), iv) the results obtained by Flowers and Itoh [T^l for the electron-electron shear viscosity in the crust region (crust, 
FI76, thick dotted curve) and v) a recent estimate by Yakovlev and Chugunov pH which accounts for the fact that in the crust 
the electrons mainly scatter off of nuclei in the crust (YC04, double dot-dashed line). The crust region is indicated by the grey 
region, and we have assumed that the core-crust transition takes place at a density of 1.7 x 10 14 g/cm 3 . All estimates (apart 
from the crust results) were determined using the simple PAL equation of state model discussed in the text. It is interesting to 
note that the Flowers and Itoh crust result dips to join our estimated rj ee as the base of the crust is approached. It is also worth 
noting the several order of magnitude discontinuity between the dominant viscosity contributions at the core-crust interface. 
Comparing the result of Yakovlev and Chugunov to the other estimates, we see that there will be a step of more than two 
orders of magnitude even when the protons are normal and provide the main scattering agent for electrons in the core. 



fraction, we also need to i) determine whether the protons and/or neutrons are superfluid, and ii) if so, quantify the 
associated suppression factors. This forces us to venture into the thorny area of nucleon superfluidity (see [2j| for a 
recent review). 

In this section, we summarise the current understanding of the relevant superfluid energy gaps and discuss the 
suppression factors we re quir e to co mple te our model. From a survey of the vast superfluid gap-literature [2(| El El 
El III El El El El El lH El El E! E3 it is clear that the determination of the relevant parameters constitute 
a serious challenge. In many ways this is not too surprising. After all, significant uncertainties concerning the 
supranuclear equation of state remain. Given the fact that the superfluid energy gaps are sensitive to the internal 
constitution of the star, we have to accept a range of possibilities at the present time. 

It is natural to adopt a strategy similar to that used in the (closely related) discussion of neutron star cooling, see 
for example the work by Kaminker and colleagues [2! 0, 143 . l43j and the recent discussion of a "minimal cooling 
paradigm" by Page et al [44j. That is, to consider the plausible range of superfluid parameters suggested by the 
theoretical studies and ask whether observations provide useful constraints on the models. Since our focus is on the 
shear viscosity in a neutron star core, the relevant observations concern pulsar glitches and potential future data for 
core fluid oscillations. In particular, it is known that gravitational waves from internal dynamics could provide an 
excellent probe of the core physics [5| • 

In principle, superfluidity in the singlet state appears if the formation of Cooper pairs from particles with opposite 
momenta and spins leads to a lowering of the ground state energy. A key parameter in any discussion of nucleon 
superfluidity is the energy gap A(fc), which corresponds to the energy needed to create a quasiparticle of momentum 
k in the superfluid E2- The gap on the Fermi surface A(kp ) can be interpreted as half the energy required to break 
the Cooper pairs. In the following we will focus on estimates of this quantity since it allows us to approximate both 
the critical temperature at which the matter becomes superfluid and the reduction factors required in our analysis of 
shear viscosity. 
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The strong interaction provides several channels in which pairing is possible. The indications are that Cooper pairs 
with zero angular momentum in a spin-singlet state should form at low momenta. However, since the effective 
interaction is expected to become repulsive near the nuclear saturation density, the corresponding energy gap will 
disappear as the density increases |3J|. As a result, because of the comparatively low proton number density one 
would expect proton x 5o pairing at densities too high for neutron pairing in that phase |27l l28j . All existing models 
suggest that the dominant proton pairing takes place well inside the core-crust transition. At the same time, the 
energetic advantage of neutron pairing in the anisotropic 3 P2 state at higher densities may be attributed to the tensor 
force H3- 

Given the practical difficulties associated with numerical many-body calculations, most studies have been based on 
the "standard" BCS theory, including more and more refined pairing interactions by adding successive corrections to 
the bare nucleon-nucleon potential |30l |. The next order beyond the bare interaction in this scheme corresponds to 
various medium effects. These are often accounted for by renormalising the single-particle energies, eg. by lowering 
the effective nucleon mass. It is by now well-known that this leads to a smaller energy gap. The key factor determining 
the effective mass is the number of interaction partners of opposite type that a given nucleon has, i.e. the number of 
neutrons per proton and vice versa. The greater the number, the greater will be the effect of the medium in reducing 
the effective mass. Hence, the effective proton mass tends to be smaller than that for neutrons, and one may expect 
the maximum gap for the proton superfluid to be smaller than that for neutrons. Similarly, a major difference 
between results for pure neutron matter and nuclear matter is due to the induced lowering of the effective neutron 
mass by the relatively small proton contaminant |29j. 

While the results for the bare interaction have converged towards a predicted maximum 1 Sq neutron gap of about 
3 MeV at k F w 0.85 fm" 1 [H HI IHHl (model A in Table 0, the higher order calculations still lead to a range 
of possibilities. Particularly important is the fact that pairing polarises the medium [2(|. After accounting for the 
polarisation effects one finds that the maximum gap has been quenched to nearly 1 MeV at kp w 0.7 — 0.8 fm" 1 
|32l l33l l37| (models a — d in Table P) . A nice explanation for this effect is provided by Lombardo and Schulze [2{| 
who show that low density polarisation effects suppress the BCS gap by roug hly a factor (4c)" 1 / 3 « 0.45. 

The Sq proton and 3 P2 neutron gaps provide further challenges. In the former case one must take into account 
the background neutron influence (eg. the polarisation terms), while in the latter one must solve the anisotropic 
gap equations (in principle ten coupled equations) in a consistent way. In the case of the proton gap the result is 
particularly sensitive to the so-called "symmetry energy" . This is expected since this parameter governs the proton 
abundance, cf. |22U3ll. The available results suggest a maximum gap for proton So pairing of approximately 1 MeV 
at k F ~ 0.4 — 0.5 fm" (models e and g in Table^l. This is roughly a factor of 3 smaller than the maximum neutron 
gap. The difference can be understood from the fact that the effective proton mass is smaller than the effective 
neutron mass. (In a pure proton medium the proton gap should, because of the charge symmetry of the nuclear 
interaction, be identical to the neutron gap in a pure neutron fluid pr|.) One would expect polarisation effects to 
further influence the proton gap, leading to a suppression (at least?) similar to that found for the 1 Sq neutron channel 
(a factor of a few) 32] . Thus the maximum proton gap may be reduced to perhaps 0.2-0.3 MeV (a level similar to 
model / in Table [IJ). 

In order to study the anisotropic 3 P2 neutron gap in detail one must extend BCS theory, and replace the single 
gap equation for the 1 Sq case by ten coupled equations. In fact, the added attraction from the tensor coupling is 
essential for the existence of superfluidity in this state (3f|. As in the case of the isotropic gaps, one finds that the 
3 P2 gap is reduced significantly (by a factor of two or so) by the lower neutron effective mass in nuclear matter. The 
available results for the 3 P2 gap illustrate the extent to which this problem remains to be understood: The calculation 
of Elgar0y et al 1351 (models k and I in TableHJ) suggests a much smaller gap than that predicted in the work by Baldo 
and colleagues |4J| (models h — j in Table QJ. The difficulties associated with superfluid pairing in the 3 P2 channel 
are exacerbated by the fact that relativistic effects come into play at the relevant densities. While the 1 S'o results 
remain largely unaffected by the inclusion of relativistic effects (mo del / in Table QJ, the associated change in the 
single-particle energies reduces the 3 P2 gap by about a factor of two [33 (model m in Table HJ). Finally, in contrast 
to the case for the 1 S'o gaps, it has been suggested that polarisation effects (which have yet to be accounted for) may 
increase the 3 P2 neutron gap. 

In order to consider variou s gap-m odels in our analysis, without complicating things excessively, we take a lead from 
Kaminker and colleagues |23l l4ll l42l l43j and represent the energy gap (at the Fermi surface) by the phenomenological 
formula 



(kF — k\) 2 (kp — k. 



V2 



A(M = A v -± e *> (19) 

(kp - kiY + k 2 (kp - k 3 j z + fc 4 

where kp is the Fermi momentum of the relevant nucleon. This expression makes sense because we know from the 
weak-coupling formula of BCS theory that A ~ k F at low densities (leaving out an exponential factor which further 
suppresses the gap as kp — ► 0). The results in the literature also indicate that the bare interaction gap function is 
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model A 


o (Mev) 


ki (far 1 ) 


k 2 (far 1 ) 


kz (fm" 1 ) 


ki (fm ] 


Reference/Comments 


A 


9.3 


0.02 


0.6 


1.55 


0.32 


[34, 36^ ^57, j59j (bare interaction) 


a 


68 


0.1 


4 


1.7 


4 


[33] 


b 


4 


0.4 


1.5 


1.65 


0.05 


[37] 


c 


22 


0.3 


0.09 


1.05 


4 


[29] (Reid potential w. m„) 


d 


2.9 


0.3 


0.017 


1.3 


0.07 


[29J (Reid potential w. polarisation) 


e 


61 





6 


1.1 


0.6 


\38j (relativistic) 


f 


55 


0.15 


4 


1.27 


4 


[271 (OPEG potential w. m, = 0.6) 


9 


2.27 


0.1 


0.07 


1.05 


0.25 


[3JJ 


h 


4.8 


1.07 


1.8 


3.2 


2 


[40] (Bonn B potential, free spectrum) 


i 


10.2 


1.09 


3 


3.45 


2.5 


[40] (Argonne Via potential, free spectrum) 


3 


2.2 


1.05 


1 


2.82 


0.6 


\40j (Paris potential, free spectrum) 


k 


0.425 


1.1 


0.5 


2.7 


0.5 


[35J (non-relativistic) 


I 


0.068 


1.28 


0.1 


2.37 


0.02 


[35J (non-relativistic) 


m 


2.9 


1.21 


0.5 


1.62 


0.5 


\38j (relativistic) 



TABLE I: Detailed parameters for our various gap models constructed from Q19I . The models are based on calculations in the 
given references and represent the current range of possibilities. Model A is for the bare interaction and is relevant in a pure 
neutron (proton) medium. Models a — d are for the 1 So neutron pairing, while models e — g correspond to the 1 S'o proton 
results and models h — m are for the 3 P2 neutron channel. 




0.5 1 1.5 0.5 



k F (fm" 1 ) k p (fm"') 

FIG. 2: The various energy gaps for the 1 So pairing, corresponding to the data in TableQ] are shown as functions of kF- The 
left panel corresponds to the neutrons, while the right panel shows the proton results. 



roughly symmetric [3 HE lU HU • This suggests that A ~ (kp — k 3 ) 2 , where the gap vanishes for kp > k 3 , for higher 
densities. The additional parameters ki and k± permit us to adapt the shape of A according to various results which 
incorporate, for example, polarisation effects. 

In Table [I] we provide data for a set of gap models based on the literature . The models have been selected 
to span the current range of uncertainty. The given parameters lead to good representations of the original results, 
although it should be noted that we have not worried too much about the numerical precision of each individual gap. 
The pairing gaps obtained from the data in Tableware shown in Figures |21 and |3J 

Having surveyed the current thinking about the nucleon superfluid energy gaps, and modelled the results in a 
useful way, we want to use the results to derive, first of all, the associated transition temperatures and then turn our 
attention to the suppression of shear viscosity. As discussed in detail by Yakovlev et al 01 ( see a ls° [30l l3lT ] for useful 
approximations), the critical temperatures can be approximated by 

f 0.5669 for 

» <^ 0.8416 for 3 P 2 (mj = 0) (20) 
1 F> { 0.4926 for 3 P 2 (\mj\ = 2) 
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FIG. 4: Critical temperatures associated with the various gap models listed in Table Q] For clarity we only show the z Pi 
neutron results for the m.j = pairing. The results for the \mj\ = 2 case would be suppressed by a factor of roughly 0.6, 
cf. I2UII . As in Figure Q the grey region represents the neutron star crust, in which there are no free protons and hence no 
corresponding T c . The inset in the right panel helps identify the two sets of gap-models we combine to describe "strong" and 
"weak" superfluidity. 



Here the effect of the anisotropy of the 3 P2 channel, for which raj represents the projection of the angular momentum 
of the Cooper pair on the quantisation axis, is apparent. (Note that, even though the \mj\ — 1 case has not been 
studied in detail it is expected to be qualitatively similar to the raj = case.) 

If we combine these approximations with our various gap models (and the data for the PAL equation of state 
from (^7J)) we obtain the results shown in Figure 2] The two panels distinguish between cases that can be refered 
to as representing "strong" and "weak" superfluidity. We also distinguish the crust and core regions. To make 
this distinction we have assumed that the transition takes place at 0.6 of the nuclear saturation density po ~ 2.8 x 
10 14 g/cm 3 , i.e. at p c « 1.7 x 10 14 g/cm 3 . It should be noted that there are no free protons in the crust, and hence 
no proton gap in that region (we are not accounting for the possibility that nuclei may exhibit pairing). It is also 
worth commenting on the fact that the models for the 1 5o neutron gap in the crust do not account for interactions 
with the crust nuclei in detail. There are many difficulties associated with this analysis which need to be addressed 
in future research. Anyway, from Figure 21 we see that superfluid neutrons are expected to be present in the crust for 
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temperatures above a few times 10 9 K. The same is true for the proton superfluid in the core. This means that, all 
but newly born (less that a month old?) neutron stars should contain both superfluid neutrons and superconducting 
protons in the 1 S'o phase. However, it is not clear whether one should expect the corresponding regions in the star to 
overlap. Figure 01 certainly gives examples of So neutron superfluids which are entirely confined to the crust. 

At this point it is relevant to mention possible constraints on the neutron star superfluidity provided by glitch data. 
By assuming that the glitches are associated with the neutron superfluid in the inner crust, one can put constraints 
on the corresponding moment of inertia required to explain the observations. This leads to a requirement that at 
least 1.5 — 2.5 % of the star's moment of inertia is in the superfluid, see for example [8j,|47j. Comparing to the models 
given in Table [IJ this constraint could mean that model d (and possibly also model e) is ruled out. Of course, this 
conclusion has many caveats and one should avoid making definitive statements given the lack of truly quantitative 
glitch models. 

As we have already discussed, the results for the 3 P 2 neutron gap are associated with a great deal of uncertainty. 
From Figure 0] we see that the critical temperature may, in fact, be as low as 10 s K. This extreme case would mean 
that the core of an accreting neutron star in a low-mass X-ray binary, for which nuclear burning in the crust is 
expected to equilibrate the temperature around a few times 10 8 K, may not contain superfluid neutrons at all. It 
should, of course, be emphasised that most of the models have significantly higher critical temperatures for the 3 P 2 
superfluidity. It would certainly be surprising, given the present results, if a neutron star with core temperature at 
the level of 10 s K did not contain a region where superconducting protons coexist with superfluid neutrons in the 3 P 2 
state. 

Next, introducing the variable r = T/T c , we learn from Ref. ^(| that the temperature dependency of the gap 
functions can be approximated by 



A . VT^T (1.456 - 0.157/ V7 + 1.764/r) for 

V= -j^jr ™{ VT=t (0.7893 + 1.188/r) for 3 P 2 [mj = 0) (21) 



Vl -t 4 (2.030 - 0.4903t 4 + 0.1727r 8 ) jr for 3 P 2 (\mj\ = 2) 



Combining these results with the critical temperatures from <|20ll we can estimate the pairing gaps at finite temperature. 
This is a key ingredient in constructing the superfluid suppression factor 1Z P1 which represents the extent to which 
the superconducting protons no longer scatter the electrons. To estimate this factor we use eq. (45) from |20| . which 
can be written 



Tip « |o.7694 + ^(0.2306) 2 + (0.07207?/) 2 + (27y 2 + 0.1476?/ 4 ) exp --^(4.273) 



\2 _1_ y2 



0.5051 



exp(4.273 - V( 4 - 273 ) 2 + V 2 ) - 1 } exp (l.l87 - ^(1-187) 2 + y 2 ) (22) 



The form of this fit to the numerical results was chosen to lead to the limits, 7Z P — > 1 as y — > and 7Z P ~ 
0.2362y exp(— y) as y — > oo. For temperatures above T c one should simply use 1Z P = 1. It is perhaps worth pointing 
out that the exact form of this suppression factor is only relevant very close to the transition temperature, and that 
an accurate model may only be important for a small subset of neutron star models. 

We need an analogous suppression factor for neutron superfluidity. In this case we have carried out a fit to the 
expressions for lZ n i and lZ n 2 given by eqn. (45) in 0- Thus we find 



Tin 



0.9543 - v / (0.04569) 2 + (0.6971y) 2 J exp (0.1148 - ^(O.IUS) 2 + 4y 2 J (23) 
obviously with y evaluated for the relevant pairing gap according to l|21|) . Given this suppression factor, we use 

r) n — > fc n r) n (24) 

for the shear viscosity due to neutron-neutron scattering. 

If we now combine all the results we have discussed, we can calculate the electron shear viscosity coefficient as a 
function of the density for various neutron star temperatures. In Figures [5J and El we show the result for strong and 
weak superfluidity, respectively. The former model combines gaps a, e and h from Table [IJ while the latter follows 
if we use gaps d, / and I. The two Figures show the shear viscosity as a function of the density for four different 
temperatures T = 10 s K, 5 x 10 8 K, 10 9 K and 5 x 10 9 K, also indicated by horisontal dotted lines in Figure 0] The 
illustrated data shows (very nicely) how the neutron viscosity becomes strongly suppressed at lower temperatures. 
One can also clearly distinguish, eg. from the transition that takes place at a density of about 7 x 10 14 g/cm in the 
last three panels of Figure [SJ the region where the protons are superconducting and the (stronger) electron-electron 
viscosity dominates. Taken together, the two figures provide a useful insight into the complexity of the neutron star 
interior and how the viscosity changes as the star matures. 
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FIG. 5: Viscosity coefficients for the "strong" superfluidity case discussed in the text. The dashed line represents the neutron- 
neutron scattering viscosity, while the solid line shows the electron viscosity. The grey region shows the extent of the neutron 
star crust. Note that, in the final panel the neutron viscosity is suppressed to levels below those shown. 




5 10 5 10 5 10 5 10 

14 3 

p (10 g/cm ) 

FIG. 6: Viscosity coefficients for the "weak" superfluidity case discussed in the text. The dashed line represents the neutron- 
neutron scattering viscosity, while the solid line shows the electron viscosity. The grey region shows the extent of the neutron 
star crust. 

These figures represent the final results of this paper. Having reached the stage where we can generate this data, 
we are now able to turn our attention to the implications for neutron star dynamics. In particular, we should be 
able to study the impact of nucleon superfluidity on neutron star oscillations in more detail than has been possible 
so far. We can also return to well-known problems like the "spin-up" problem and re-assess the viscosity due to the 
core-crust interaction and the associated Ekman layer. 
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IV. FUTURE CHALLENGES 

Hoping to stimulate further work on the many relevant issues, we have discussed the effects of superfluidity on 
the shear viscosity in a neutron star core. As should be clear from our analysis, an understanding of this problem 
requires input from several different areas of research, most notably nuclear physics and relativistic astrophysics. 
Although the results we have collated may not be original, we are not aware of any previous study which combines 
the data into a consistent description of the shear viscosity in a superffuid neutron/superconducting proton mixture. 
This is in sharp contrast to the many discussions of the role of superfluidity in neutron star cooling, see for example 
[23L liH I42I 143 . l44j . We believe that this paper provides a useful update on the core viscosities which should be 
valuable for future modelling of neutron star oscillations and associated instabilities [3( ■ 

We have provided a simple model for the electron viscosity which is relevant both when the protons form a normal 
fluid and when they become superconducting. This result explains (as in why it is natural that proton supercon- 
ductivity leads to a significant strengthening of the shear viscosity. It also clarifies some confusion associated with the 
description of Cutler and Lindblom |14j . In particular, it is clear that the superfluidity of the neutrons is not the key 
factor which leads to electron-electron scattering becoming the main shear viscosity agent. Rather, it is the fact that 
the onset of superconductivity suppresses the electron-proton scattering. This is an important point. As can be seen 
from Figure ^ the electron-electron shear viscosity is not too different from the result for neutrons scattering off of 
each other. This means that, in the temperature range where shear viscosity dominates, the damping of neutron star 
oscillations will be quite similar (modulo multi-fluid effects) in the extreme cases when i) the neutrons and protons 
are both normal, and ii) when the neutrons and protons are superfluid/superconducting, respectively. The contrast 
with the case when the neutrons are superfluid and the protons normal (and viscosity is dominated by rj ep ) is clear 
from Figure^ This is an interesting observation because it shows that proton superconductivity (or rather absence 
thereof) could have a significant effect on the dynamics of a neutron star core. 

Our results are in a form which permits the use of data for any given modern equation of state. In contrast to the 
often used expression obtained by Cutler and Lindblom |14| , our final formulas are explicitly dependent on the proton 
fraction (as well as the total mass density). This is a crucial difference because of the simple fact that the result in 
was obtained for the now seriously outdated equation of state derived by Baym, Bethe and Pethick 0] m 1971. 
We can obtain an indication of the likely "error bars" associated with the electron shear viscosity from l|ltj|) which 
suggests that r) e ~ x p ' . As the proton fraction varies by up to perhaps a factor of two for different equations of state 
we can expect the viscosity coefficient to be uncertain at the level of (at least) a factor of three. 

Our discussion also highlights a number of challenges that must be met if we are to make further progress. First of 
all, it is clear that an understanding of the various superfluid gaps is crucial. We believe that the results summarised in 
Table[U i.e. the "phenomcnological" models we have constructed from a range of results in the relevant literature (in 
line with the philosophy adopted in [2^, 0,E3>E3l)j provide a useful survey of the theoretical range of possibilities. 
These models are also readily included in quantitative studies like that discussed in this paper. Of course, many 
uncertainties remain concerning the superfluid parameters. Most important are the suppression factors for the various 
scattering rates which are needed to evaluate the viscosity coefficients (1Z P and lZ n in our analysis). Any detailed 
discussion of these coefficients is certainly valuable. Especially since it allows us to make the analysis of temperatures 
near the critical temperature T c quantitative rather than qualitative (as in the present discussion). 

In addition to these problems, there are many related challenges. Although it is commonly acknowledged that 
superfluid components play a crucial role in neutron star dynamics (eg. in pulsar glitches |2|, [7J, |8| ) , our understanding 
of the multi-fluid aspects and the relevant dissipation mechanisms can be improved considerably. That this area 
provides exciting possibilities is nicely illustrated by the demonstration that a two-stream instability may operate 
in superfluid systems with relative flow |48f. Of particular importance for the physics in the neutron star core will 
be the superfluid vortices (and magnetic fluxtubes in the case of a superconductor), leading to i) dissipation due to 
mutual friction l49l| and ii) potential vortex pinning to the crust (see, for example, the discussion in |5(j). As 
the parameters required to model dynamical superfluids (like the entrainment coefficients |l2t \5l\) are not yet well 
constrained, this is a research area in its infancy. 

Yet another illustration of this fact is provided by the recent discussion following the likely observation of free 
precession in PSR B1828-11 |3,|52|. Of particular relevance for the present paper is Link's suggestion [l(| that the 
interaction between neutron vortices and proton fluxtubes in a type-II superconductor would lead to a strong coupling 
between the two fluids and prevent long-period precession. As an alternative, Link suggests that the protons could 
form a type-I superconductor, from which the magnetic flux would be expelled. This idea prompted a recent analysis 
by Buckley et al |53| | , who find that the interaction between neutron and proton Cooper pairs may indeed lead to the 
neutron star core becoming a type-I superconductor. Exactly how this would affect the various viscosity estimates, 
and the dynamics of the neutron star core, is unclear at this point. 

Models of the neutron star crust almost exclusively assume that it can be described as a regular bec lattice. At 
some level this will be an oversimplification. In particular, it is likely energetically favourable for the nuclei to have a 
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variety of "shapes", ranging from droplets to rods to plates, as the core-crust phase-transition is approached |54|. An 
effective description of these exotic phases may differ significantly from the standard picture. In addition to this, one 
should account for the presence of superfluid neutrons throughout the bulk of the crust. The relevant parameters for 
this superfluid phase must be considered as poorly known since the various gap-calculations do not (yet) account for 
the presence of the nuclei. This constitutes a difficult theory problem, but the result could be crucial for a number 
of scenarios ranging from the pulsar glitches to gravitational-wave driven mode-instabilities. An indication of how 
different the parameters of the crust superfluidity may be compared to those in the core fluid is provided by the work 
of Carter, Chamel and Haensel |55l 1561 157| . Their analysis of the entrainment near the base of the crust shows that 
the results differ greatly from those typically used to model the outer core. The essential difference is in the effective 
mass of the neutrons. As in the other cases, the dynamical repercussions of this result remain to be investigated. 

Finally, let us turn to the physics of the neutron star inner core. It is generally expected that, at densities above 
a few times nuclear saturation density, it will be energetically favourable for exotic phases of matter to be present. 
The two most commonly considered cases concern hyperons and deconfined quarks. In the first of these cases it has 
been pointed out that hyperons provide an efficient refrigerant, since direct URCA reactions will be in operation. 
These reactions would cool the star much faster than indicated by the observational data. The only way to avoid 
contradictions with the observations is to appeal to hyperon superfluidity |58| , which would quench the reactions and 
slow down the cooling rate. It is also worth pointing out that the fast nuclear reactions in a hyperon core are thought 
to lead to rapid damping of oscillations due to the resultant bulk viscosity [5^. l60| . This effect would be significantly 
suppressed by superfluidity 61]. In the context of the present paper it is relevant to make two comments. First of 
all, it is known that an increase of the hyperon population tends to lead to a significant depletion of the electron 
number density. If we suppose that all the other components (neutrons, protons and hyperons) are superfluid, then the 
electrons may provide the only channel for shear viscosity. If their number density is depressed to say x e ~ 10 -3 then 
we would estimate that the shear viscosity is something like five orders of magnitude weaker than that illustrated in 
Figure^ The dynamical effects of this could be important. Secondly, it is worth pointing out that neither the many 
possible vortex- vortex and electron- vortex interactions, nor the multi-fluid aspects of superfluid hyperon mixtures 
have yet been addressed. 

The situation is quite similar for the deconfined quark core. The mixture of s, d, and u quarks is thought to be 
able to exhibit a variety of "colour superconducting" pairings |62j . Whether this leads to multi-fluid dynamics is not 
clear, although it seems likely for at least some of the possible phases. There have not yet been many attempts to 
estimate the relevant transport coefficients. Yet, it is clear that such studies are of great importance. Especially since 
they may help unveil the true ground state of matter. In this context, it is worth mentioning an interesting result. 
Madsen [63| has combined estimated viscosity coefficients for strange stars with the radiation reaction results for the 
r-mode instability to demonstrate that, if the fastest spinning pulsar were a pure colour-flavour locked quark star then 
it ought to be unstable (see [64] for an alternative). One interpretation of this result is that this particular form of 
quark pairing is not present in the bulk of a compact star. There are of course many caveats to this statement, but 
it is nevertheless relevant. In particular, it demonstrates the need to put current and future models of exotic neutron 
star physics in an astrophysical context and ask whether observations provide useful constraints on the theory. It is 
particularly important to consider this at the present time when gravitational-wave observations of compact object 
dynamics seem a definite possibility. 

Taking the situation at face value, neutron star (astro) physics is a vibrant area of research which provides a number 
of exciting challenges for the future. In order to meet these challenges we need a closer dialogue between researchers 
in different fields, like nuclear physics and relativistic astrophysics. We hope that this paper will help stimulate such 
discussion. 
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